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Abstract. The stability properties of a class of dissipative quantum mechan- 
ical systems are investigated. The nonlinear stability and asymptotic stability 
of stationary states (with zero and nonzero dissipation respectively) is investi- 
gated by Liapunov's direct method. The results are demonstrated by numerical 
calculations on the example of the damped harmonic oscillator. 



1. Introduction 

The study of dissipative quantum mechanical systems is of widespread interest 
in different areas of physics HE]. In several approaches the microscopic systems 
are embedded in an environment and the dissipation is given with a detailed de- 
scription of that assumed macroscopic background. However, the structure of the 
environment is principally unobservable or unknown, moreover, its details are not 
very important but only some basic properties. Several developments of the sto- 
chastic approach to decoherence and quantum dissipation can be interpreted as a 
search after these important properties and get rid of the details of the background. 
The unacceptable negative densities in the soluble Caldeira-Leggett model |3] can 
be cured by additional damping terms ^ according to the requirement of positivity 
of the irreversible dynamics ,5,. However, complete positivity alone does not give 
automatically physically meaningful results, there are indications that some further 
refinements are necessary [HJ [7| . 

On the other hand, the approaches from the phenomenological, thermodynamic 
side |H1 El lini ITT] are fighting with the right interpretation and properties of ther- 
modynamic quantities in quantum systems. 

The situation is further complicated with the fact that there is no common agree- 
ment in the nature of quantum dissipation, there is no unique quantum version of 
the simplest dissipative classical systems like the damped harmonic oscillator. The 
problem is that the lack of a uniformly accepted variational principle prevents to 
use normal canonical quantization |12| . Hence quantization can be accomplished in 
several different ways and the quantizations based on different variational princi- 
ples are not equivalent at all, moreover show peculiar physical properties. The best 
example and a kind of parade ground of the different methods is the mentioned 
quantum harmonic oscillator with the simplest possible damping I14j . 
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In this paper we approach dissipation in quantum mechanics from a different 
point of view and, instead of the von Neumann equation, we investigate the prop- 
erties of some simple generalizations of the Schrodinger-Langevin equation, where 
the damping force is negatively proportional to the velocity 

at 2m 2 ip* 

The above equation is a quantized form of Newton equation of a classical mass- 
point with mass m, moving in a potential U and with damping term proportional 
to velocity 

(2) mx = rTiv = — VC/ — fcv. 

One can give several derivations of Schrodinger-Langevin equation based on dif- 
ferent assumptions ^lEJEIElElElEni and the equation has several appeahng 
properties. It can be derived by canonical quantization, the dissipation term has 
a clear physical meaning, the equation preserves the uncertainty relations, and all 
the stationary states of a quantized Hamiltonian system are necessarily stationary 
states of the corresponding dissipative quantum system. The later property can 
be true in more general systems, too [2Ij. At the first glance this last fact seems 
to be in contradiction with the apparent instability of excited states |22| and in- 
cited some criticism therefore it requires a more detailed investigation that can 
be expected partially from stability investigations. The relationship between the 
Schrodinger-Langevin equation and the master equation approaches is important 
also from fundamental theoretical point of view, as it was indicated for example in 

El 

There are some known exact solutions of equation Q for the cases of free motion, 
motion in a uniform field, and for a harmonic oscillator in one dimension [251 12(ij . 
These exact solutions have the remarkable property that the stationary states of 
the system are asymptotically stable. 

We will show that the above properties are valid not only for these known exact 
solutions of the damped harmonic oscillator but in general for any kind of initial con- 
ditions, moreover not only for the oscillator, but in case of (almost) arbitrary poten- 
tials. The stability investigation is based on the hydrodynamic model of quantum 
mechanics, which recently gained a renewed interest because of the development of 
powerful numerical codes of computational fluid dynamics j27[ I28L I29L 1301 1ST] , and 
of its relationship to some quantum field theories as non-Abelian fluid dynamics 
|H2j , in perturbational cosmological calculations curing some instabilities of the Eu- 
ler equation PSI, etc. In the last section the above properties will be demonstrated 
numerically on the traditional parade ground of dissipative quantum mechanics, on 
the example of the damped harmonic oscillator. 



2. Stability properties of the damped Schrodinger equation 

Let us consider a classical mass-point with mass m, moving in a potential U , 
governed by the following Newton equation: 

d^x dv „ 

(3) m^=..- = -VC/ + F„ 

where U is the potential of the conservative part, Fc — —VU, of the total force. 
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It is well-known that the equation of motion of the corresponding quantized 
system can be given when Fd = and is the Schrodinger equation, with k = 
0. The Schrodinger equation can be transformed into various thought provoking 
classical forms via the Madelung transformation 

and by introducing the so-called hydrodynamic variables, defining the probability 
density p and a velocity field v as 

(4) p := i?2 = 

(5) V :== — VS* = -— Vln-7-. 

m 2m ip* 

Here the star denotes the complex conjugate. With these variables the real and 
imaginary parts of the Schrodinger equation give 

(6) g+v.(pv) = 0, 

(7) m— = -V (my + + , 

with the condition that the vorticity of the "fluid" is zero: 

(8) cj = Vxv = 0. 
Here 

Ua = = Ap - ^ ' ' 

^ 2m R Amp \ 2p 

is sometimes called the quantum potential, and the equations (0-0 ^-^'^ governing 
equations of the fields (p, v). In particular, © is the balance of "mass" or proba- 
bility density. The second equation Q is the Newtonian equation for a point mass 
moving in superposed normal and quantum potentials. The governing equations 
are more apparent if we transform them into a 'comoving' form: 

(9) p + pV-v = 0, 

(10) mv = -V [Uq -f U) . 

Here the dot derivative denotes the substantial time derivative, 9t -f v ■ V. It 
is easy to see that the above equations preserve the probability and vorticity. In 
the hydrodynamic formalism it is apparent that any other kind of dissipative forces 
that would destroy the condition of zero vorticity would destroy the connection 
to the single Schrodinger equation and initiate coupling to the electromagnetic 
field or/and the development from pure into mixed states. It is straightforward to 
transform the above system into a true hydrodynamic form, because the quantum 
potential can be pressurised that is, the corresponding force density can be written 
as a divergence of a quantum pressure tensor Pq-. 

pVUq^y-Pq, 

where the quantum pressure is not determined uniquely from this condition. It can 
be chosen as 

ft2 / VpoV(o' 
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or as 



p = 



(V o Vjp — 



-(VoV)lnp, 



P 

or in some other form that differs only in a rotation of a vector field. Here o is 
the traditional notation of tensorial product in hydrodynamics, and / is the second 
order unit tensor. Let us remark that thermodynamic considerations result in a 
unique quantum pressure j34| . Therefore, the second equation of the above system, 
(jlOjl , can be written in a true Cauchy form momentum balance as well, introducing 
the mass density p — mp, 

(11) g + v(/5v)^0, 

(12) ■^^ + V(pvov + P,) =-pVC/. 

All these kinds of transformations of the Schrodinger equation are well- known 
in the literature with different interpretations. The Newtonian form was investi- 
gated and popularized e.g. by Bohm and Holland jSSlESI- The hydrodynamic form 
was first recognized by Madelung |37| ) and extensively investigated by several au- 
thors, e.g. ,38».^40|- Moreover there is also a microscopic-stochastic background 
|41[l42lllSll44| . giving a reasonable explanation of the continuum equations. Inde- 
pendently of the interpretation, one can exploit the advantages of the hydrodynamic 
formulation to solve the equations or investigating its properties. 

First of all let us remark that the quantum mechanical equation (|10|l can be 
supplemented by forces that are not conservative in classical mechanics but still 
admit a hydrodynamic formulation: 

(13) mv = -V ([/,-(-£/)+ Fd(v). 

The corresponding balance of momentum is 

(9pv) 
dt 

It is easy to check that all force fields of the form Fd = ~k{S)v are rotation free, 
therefore do not violate the vorticity conservation and the corresponding quantized 
form can be calculated by canonical quantization with the help of the velocity 
potential S The hydrodynamic transformation gives the quantized form of the 
damping force expressed with the wave function as Fd = _l!^M£l\7 in Moreover, 
the Schrodinger equation is transformed into the Schrodinger-Langevin form 
Let us mention here that the nonlinearity of equation means only a practical, 
not a fundamental problem '2U'. 

The real stationary solutions of the above system are those where the substantial 
time derivatives are zero. Thus we can eliminate the virtual effect of the motion of 
the continuum, and the stationary solution will not depend on any kind of external 
observers, stationarity becomes a frame independent property. Then, the nontrivial 
{p ^ 0) stationary solutions (ps,Vs) coincide with the stationary solutions of the 
Schrodinger equation without damping, and in the hydrodynamic language are 
characterized by the conditions, 

(15) V,, = 0, 

(16) Uq{ps) + U ^Es^ const. 



(14) ^^-HV(pvov + P,)-p(-V{/ + Fd). 
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For the sake of simphcity, in what foUows we choose units in which h — 1 and 
m — 1 and which make all quantities dimensionless. 

In the following we investigate the stability of the stationary solutions (|15ll - llt)|l 
of the equations Hll() and (|14|l . 

Let us assume that the functions (p, v) are two times continuously differentiable 
and the density ps exponentially tend to zero as |x| goes to infinity 



(17) 



lim ps(x)el''l = 0. 

|x| — *oo 



Furthermore the dissipative force obeys the following inequality for all possible 
solutions of the above equations 1)11(1 and ((14|l 



(18) 



pv-FddV < 0. 



Here equality is valid only if p-Vg ~ 0. 
With these conditions we show that 



(19) 



2p 



dV 



is a Liapunov functional for the stationary solution in question. 

First we show that, (|19|l is decreasing along the solutions of the differential 
equation. Indeed, the first derivative (variation) of L for {pp (5v) is 



5L: 



/9V(5v 



U^Uq-EAbp 



dV 



Vlnp 



(5/9 ■ dA, 



where the surface for the second integral is a sphere with the radius increasing 
to infinity. That term zero at the corresponding equilibrium (stationary) solution 
under the asymptotic condition ((17|l . Therefore the derivative according to the 
differential equation is 
(20) 



dL 
dt 



pw 



Ua + U-E, 



Vlnp V-(pv) 



• dA + / pv-FddV <Q, 



because the first surface integral vanishes according to the asymptotic conditions, 
and the second integral is negative due to (fl^ . 

Next, the definiteness of the second derivative of ((19|l is investigated. Evaluating 
the second variation between (^Sp (5v) and (^S'p S'v) that yields 



(21) / {Sv 6p VSp) 




(<5'v 6'p Vd'p) dV. 



4p2 4p 

We can see that the minors of the above matrix under the integration are 



(22) 



A3^-^. 

4p 



The matrix in H21|l is positive semidefinite in case of stationary solutions, there- 
fore the linearized system is marginally stable. 
Remarks: 
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- The continuity and asymptotic conditions are valid for stationary solutions 
in case of several classical potentials (harmonic, Coulomb, etc.). More- 
over, one can find different (weaker) asymptotic conditions to ensure the 
differentiability of L. 

- The physical meaning of the Liapunov function (|f 9|l is that it is essentially 
the expectation value of the energy difference of the given and the stationary 
states. The decreasing property of L is simply the decreasing of energy in 
our dissipative system. The Liapunov function can be expressed with the 
wave function as 

(23) L(^) = - 1 (^(V^)' + ^ W*)' + {U- E,)^l^r) dV 

- According to the structure of the second derivative of L, given in (|^ . 
the function space of (p, v) functions looks like a Soboljev space with our 
asymptotic conditions and norm 

ii(p,v)iis=|y' [p'+^^+ivpf] d^y^' 

- If Fd = then the equilibrium solution is stable with the above condi- 
tions. We can see that the dissipative, damping force makes the station- 
ary solutions attractive, therefore we can expect asymptotic stability (the 
conditions can be rather involved see e.g. |45j). In ordinary quantum 
mechanics without damping forces the stationary solutions can be only sta- 
ble. Furthermore, the above condition can serve as a definition of damping 
(dissipative) forces, i.e. we can say that a force is damping (dissipative) if 
pv • Fd < and / pv • Fd dV < for any configuration p, v. 

- The marginal stability indicates a rich structure of possible instabilities and 
therefore the need of further stability investigations. On the other hand 
a rich stability structure is expected considering the infinite number of 
stationary states. Moreover, in more than one dimension and if the domain 
of the functions (p, v) is not the whole space (e.g. in case of a Coulomb 
potential), the perturbations destroying the rotation free streams requires 
to consider Casimir functions for the stability investigations |45l 1^ . 

- The given conditions are capable to estimate the basin of attraction of the 
different quantum states. The physical role of dissipation is to relax the 
process to the lowest energy element within a basin of attraction. Moreover, 
the above local (!) statement regarding the formal stability of stationary 
states does not exclude the possibility that for a given initial condition the 
whole set of stationary solution is the attractor of the quantum dynamics. 
On the other hand transition from one basin to another possibly can be 
induced by some additional interaction introduced in the system. 



3. Stability properties of the harmonic oscillator 

As a demonstration of the previous stability properties we investigate the one 
dimensional dissipative harmonic oscillator, with constant damping coefficient. In 



STABILITY OF STATIONARY SOLUTIONS OF THE SCHRODINGER-LANGEVIN EQUATION7 



this case the hydrodynamic equations are simphfied into the foUowing form 
(24) 

(25) 



dp 


d{pv) 


at 


dx 




dv 


d 


m — 
2 


dx 




kv. 



Let the initial condition be a standing wave packet (winit = 0), having the shape 
Pinit — P(o,o.o.5)j ^ind with its maximum shifted from the center of the osciUator 
potential to the left by two units. We will investigate the time development of this 
initial distribution in a potential with elasticity coefficient D = 0.02 and different 
damping coefficients. In the numerical calculations the equations have been solved 
with a modified leapfrog method, adapted to the presence of the nonlinear hydro- 
dynamic term and for dissipative forces. The asymptotic conditions at the infinity 
were considered by extrapolating the internal velocity values at the boundaries. 

On the first two figures one can see the time development of the probability 
distribution and the current, respectively. In this calculations the damping was 
large, with a coefficient k = 1. On figure^the center of the oscillator is at point 41 
and the solutions starting from P(o,o.05) are given in every 5 time units from to 
100. One can observe that the distribution tends toward the stationary solution, 
denoted by dots. At the same time, the probability current goes to zero as one can 
see on figure El It is interesting to observe that, in the beginning, there is a space 
interval where the velocity is negative and later the direction of the motion become 
everywhere positive. 




Figure 1. k = 1, large damping. Probability distribution of the 
position, < = 0, 100, in every 5 time units. 

At the calculations leading to the next two figures the damping was smaller, 
with a coefficient k = 0.1. On figure|3|the center of the oscillator is at point 21 and 
the solutions are given in every 1 time units from to 20. One can observe that 
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Figure 2. fc = 1, large damping. Probability current of the posi- 
tion, t = 0, 100, in every 5 time units. 



the trend to equilibrium is dilferent from that of the previous case. A new period 
of a damped oscillation starts, as one can conclude from figure 0] of the probability 
current. At figures ISlHI there are the undamped solutions with fc = 0. The center 

P 



0. 25 
0. 2 

0. 15 
0. 1 

0. 05 




Figure 3. fc = 0.1, small damping. Probability distribution of 
the position, t = Q, 20, in every 1 time units. 

of the oscillator is at point 21 and the solutions starting from jO(o,o.05) ^-re given in 
every 0.95 time units from to 19. 



STABILITY OF STATIONARY SOLUTIONS OF THE SCHRODINGER-LANGEVIN EQUATION) 



0.03 - 



0.02 - 



0.01 - 



12 . 5 1 




27.5 30 



-O.OlL 



Figure 4. k = 0.1, small damping. Probability current of the 
position, t — 0, 26, in every 1.3 time units. 




17. 5 



Figure 5. fc = 0, zero damping. Probability distribution of the 
position, t ~ 0, 19, in every 0.95 time units. 



On figure one can see the distance from the stationary solution in the 
three cases as a function of time. Remarkably, the starting oscillations are clearly 
nonlinear. The new, starting periods indicate clearly that the norm is not the 
suitable one for stability considerations. The expected asymptotic stability requires 
a monotonous tending to the stationary solution. 
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Figure 6. fc = 0, zero damping. Probability current of the posi- 
tion, f = 0, 19, in every 0.95 time units. 

Dist 



0.035 r 




20 40 60 80 100 



Figure 7. L? distance from the equihbrium. k = 0,0.1, 1 at the 
dotted, continuous and dashed Hnes respectively. 



4. Conclusions 

The stabihty properties of the stationary sohitions of the Schrodinger equation 
were investigated. It has been shown that the stabihty of stationary states can be 
treated by Liapunov's direct method, where the expectation value of the energy is 
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a good starting point for stability investigations. However, the semidefiniteness of 
the second derivative at the stationary states and the infinite number of stationary 
solutions of the equation indicates a rich stability structure. These properties are 
more apparent in the chosen hydrodynamic model of quantum mechanics, enabling 
the use of the methods and a direct comparison of the well understood hydrody- 
namic stability problems. Moreover, the hydrodynamic analogy shows that in one 
space dimension we do not need Casimirs therefore the energy is also sufficient for 
exact and physically relevant results. 

We have shown that in the nondissipative case, for the Shrodinger equation one 
can expect only stability. Moreover, the above properties were demonstrated for the 
damped harmonic oscillator. This property indicates that an arbitrarily perturbed 
nondissipative quantum system cannot move from a given stationary state to an 
other one, it will oscillate around the perturbed state. 

Finally we would like to emphasize again that the stability results are indepen- 
dent on any interpretation and could be formulated with the help of wave functions 
or stochastic processes, too. However, in case of dissipation the hydrodynamic for- 
mulation has the advantage of a clear nonequilibrium thermodynamic background 
with well established stability results and clean concepts of dissipation. This kind 
of distinction of dissipation types can be crucial if one would like to compare and 
extend the above calculations to successful and experimentally tested frictional 
models as e.g. the one of Gross and Kalinowski for heavy ion collisions |47l I48| . 
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